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Abstract 

The effective theory for the dynamics of hot non-Abelian gauge fields with spatial 
momenta of order of the magnetic screening scale g 2 T is described by a Boltzmann 
equation. The dynamical content of this theory is explored. There are three relevant 
frequency scales, gT , g 2 T and g 4 T, associated with plasmon oscillations, multipole 
fluctuations of the charged particle distribution, and with the non-perturbative gauge 
field dynamics, respectively. The frequency scale gT is integrated out. The result is 
a local Langevin-type equation. It is valid to leading order in g and to all orders in 
log(l/g), and it does not suffer from the hard thermal loop divergences of classical 
thermal Yang-Mills theory. We then derive the corresponding Fokker-Planck equation, 
which is shown to generate an equilibrium distribution corresponding to 3-dimensional 
Yang-Mills theory plus a Gaussian free field. 
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At sufficiently high temperature T the running gauge coupling g = g(T) in non- 
abelian gauge theories is small. Nevertheless, long distance modes of hot non-abelian 
gauge fields are strongly coupled. This leads to the phenomenon of magnetic screening 
on the length scale (g 2 T) _1 . 2 It distinguishes non-abelian from abelian plasmas in 
which magnetic fields can be correlated over arbitrarily large distances. 

The strong coupling is caused by the large amplitudes of the infrared gauge field 
modes. It can be easily understood as follows. The long wavelength (A 3> T _1 ) 
modes obey the classical Rayleigh- Jeans law which states that the energy density is 
proportional to T. Then, by dimensional analysis, the magnetic energy density B 2 due 
to wavelengths of order A must be of order TA~ 3 . This corresponds to ~ T" 1 / 2 ^- 3 / 2 
and to a vector potential ~ T l l 2 \~ l l 2 . The Yang-Mills equations of motion are 
non-linear since they contain covariant derivatives. When A approaches ((? 2 T) _1 , both 
terms in the spatial covariant derivative D = V — gA become of order g 2 T, and the 
term gA can no longer be treated as a perturbation. 

This breakdown of perturbation theory was first noted in the thermodynamics of hot 
non-abelian gauge fields [1]. The non-perturbative physics can be ascribed to euclidean 
pure Yang-Mills theory in 3 dimensions which is obtained from the 4-dimensional ther- 
mal field theory in the imaginary time formalism by integrating out all modes with 
characteristic wavelengths shorter than (g 2 T) -1 [2]. The 3-dimensional theory can be 
easily treated non-perturbatively on a lattice. 

Dynamical quantities, which are determined by the real (Minkowski) time evolution 
of the A ~ (g 2 T) _1 modes, are more difficult to deal with. An important example is 
the so called hot sphaleron rate. In the standard electroweak theory it determines the 
rate for anomalous baryon number violation and is therefore a crucial ingredient in 
particle physics models which try to explain the baryon asymmetry of the universe [3] . 

Fortunately, even for dynamical quantities one can use perturbation theory to inte- 
grate out short distance modes to obtain an effective theory for the non-perturbative 
long-distance dynamics. The first step is to integrate out "hard" physics associated 
with virtual momenta of order T. At leading order in g one obtains the so called hard 
thermal loop effective theory [4, 5]. It describes gauge fields with A ^{gT)' 1 interacting 
with classical colored particles. These particles correspond to quanta of the |k| ~ T 
field modes which have virtual momenta of order gT or less. In a second step one 
can integrate out the degrees of freedom associated with /c , |k| ~ gT. The resulting 
effective theory is described by the classical field equations of motion [6] 

D,F^ = mlW\ (La) 
(C + vD)W = v-E + £, (l.b) 

where F£ v = d^A u a — d h ' + g] a } )C A^A v c is the non-abelian field strength tensor, and 
t> M = (l,v). The field W(x, v) represents the fluctuations of adjoint color charge due 
to hard particles with 3-velocity v, v 2 = 1 3 . For other approaches leading to the 

2 The units are chosen such that % = c = fee = 1- Furthermore, it is assumed that fermion masses 
and chemical potentials can be neglected relative to T. 

3 For notational simplicity we write all expressions in 3 spatial dimensions. One has to keep in 
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Boltzmann equation (l.b) which do not make use of the hard thermal loop effective 
theory, see [7]- [10]. The current on the rhs of Eq. (l.a) is given by 

W(x) = J v"W(x,\), (2) 

where J v = / dfl v /(4:7i), times the square of the leading order Debye mass inn oc gT. 
The lhs of Eq. (l.b) contains the linear collision term 

CW(x,v)= [ C(v,V)W(x,V), (3) 

which is due to scattering of hard particles with velocities v and v' with a momentum 
transfer of order gT. The fields in Eq. (1) describe fluctuations close to thermal equilib- 
rium on length scales larger than (gT) -1 . The collision term, breaking time reflection 
invariance, is dissipative. It is accompanied by the Gaussian white noise £, which has 
vanishing expectation value. Its only non-trivial correlation function is 

2T 

{U*,v)ZbWy)) = —C(v,V)6 ab 6 4 (x - x 1 ). (4) 
m D 

The collision kernel C(v,v') is of order g 2 T. Its precise form is not important for 
the purpose of this letter. We will only use the fact that C commutes with rotations 
in v-space. Thus if W(x, v) is expanded in spherical harmonics, the collision operator 
becomes diagonal and its eigenvalues q only depend on I, i.e., 



J C(v,V)^W lm (x)Y lrn (v') = '£c l W lm (x)Y lrn (v). 
Im Im 



(5) 



C consists of a piece which depends logarithmically on the infrared cutoff for the 
|k| ~ gT modes and of a cutoff independent part. The / = eigenvalue c vanishes. 
For the logarithmic part this was shown in Ref. [6]. In Ref. [11] it was found that 
this holds for the complete cq if one uses dimensional regularization to define C. The 
vanishing of Cq ensures that Eq. (l.b) is consistent with Eq. (l.a), i.e., that the current 
on the rhs of Eq. (l.a) is conserved. The complete / = 1 eigenvalue c\ in dimensional 
regularization was explicitly calculated in Ref. [12]. 

At leading order in log(l/g) _1 Eq. (1) can be approximated by the Langevin equation 

[6] 

D x B = 7 E + C, (6) 

where E l = F l ° and B l = —^^F^ are the non-abelian electric and magnetic fields. 
The damping coefficient or color conductivity 7 is proportional to T/log(l/g). The 
Gaussian white noise C, satisfies 

(Q(x)(i(x')) = 2 1 T6 ab 6 i i6(x-x'). (7) 



mind that the effective theories discussed here, except Eq. (6), require regularization, for example by 
a continuation to d = 3 — 2e dimensions. 
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Therefore it keeps the gauge fields in thermal equilibrium at temperature T. Eq. (6) 
implies that the characteristic frequency of the gauge fields is of order log(l/g)g 4 T [13]. 

One motivation for integrating out the momentum scale |k| ~ gT was that the Hard 
Thermal Loop effective theory contains UV divergences which cannot be removed by 
renormalization. Thus simulations of the non-perturbative gauge field dynamics using 
the HTL effective theory do not have a continuum limit. This problem still persists in 
(1), however. If one would try to take the continuum limit of (1) by sending the UV 
cutoff to infinity one would encounter precisely the same divergences, simply because 
the collision term and the noise, which distinguish Eq. (1) from the HTL effective 
theory, can be neglected at very large momenta as they do not grow as fast as the 
derivative terms. 

The effective theory (6), on the other hand, is UV finite [7]. Thus it is well suited for 
non-perturbative lattice simulations. It was used to compute the hot sphaleron rate 
by Moore [14]. In Ref. [15] it was extended to account for a Higgs field with a thermal 
mass of order g 2 T. 

Recently Arnold obtained a non-local Langevin equation for the k ~ g 4 T dynam- 
ics which is valid to leading order in g and to all orders in log(l/(?) _1 [16]. Arnold 
and Yaffe showed that it can be used to systematically improve the theory (6) in a 
perturbative expansion in log(l/0) _1 [12]. They found that Eq. (6) is still valid at 
next-to-leading order in log(l/g) -1 if one includes a next-to-leading log correction in 
the color conductivity 7 [12] . With this correction Moore's result for the hot sphaleron 
rate [14] agrees surprisingly well with different simulations of the hard thermal loop 
effective theory [17]- [19], which include all orders in log(l/(?) _1 but which do not have 
a continuum limit [20]. 

The purpose of this letter is to fully explore the dynamical content of Eq. (1). It 
will be shown that, in addition to the well known plasmon oscillations and the non- 
perturbative gauge field dynamics, Eq. (1) describes fluctuations of multipole moments 
of W with a characteristic frequency of order g 2 T. Then we will obtain a generalization 
of Eq. (6) by integrating out the physics of plasmon oscillation which is characterized 
by the frequency scale gT. The result, Eq. (16), is a Langevin equation which is local 
in space and time. In contrast to Eq. (6) and to the Langevin equation of Ref. [16] it 
contains two different frequency scales, g 2 T and g A T. 

In the following power counting estimates logarithms of g will be ignored. All ap- 
proximations will be valid at leading order in g and all orders in \og(l / g)~ l . Sometimes 
it will be convenient to write the scalar, vector, and multipole components of W sepa- 
rately. The latter will be denoted by W, so that 

W(v) = IU° + 3v- W + W(y), (8) 

where the factor 3 simply follows from the definition (2) and from / v t)V = 

The only length scale in Eq. (1) is set by C' 1 and the magnetic screening length 
which are both of order ((? 2 T)~ 1 . We have already seen that the covariant derivative 
D is then of the same order of magnitude as the ordinary derivative. Therefore one 
can obtain important information about the dynamical content of Eq. (1) already by 
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considering the linearized equations of motion. We Fourier transform them, 

ik°E + ikxB = m^W, (9.a) 
ik-E = mlW°, (9.b) 
(C-ivk)W = v-E + e, (9.c) 

and consider only |k| ~ g 2 T. The magnetic field can be eliminated from Eq. (9. a) using 
k x E = A;°B, 

ik E + -^k x k x E = m 2 D W. (10) 

First consider the case ko ^> |k|. Then one can neglect the second term on the lhs 
of Eq. (10), so that 

E=-^W. (11) 

Combining Eqs. (9.b) and (11) one finds that W° ~ k ■ W//c° <C |W|. Thus one can 
neglect W° in Eq. (9.c). Since k ^> C ~ g 2 T one can neglect the collision term and 
the noise in Eq. (9.c) and one can approximate v ■ k ~ ko, which gives 

2 

k (3v • W + W) = • W. (12) 

If one multiplies Eq. (12) with v and integrates over v, W drops out and one obtains 

(kl - \ml) W = 0. (13) 

Thus W oscillates with the plasmon frequency u p \ = 4^mD- The plasmon oscillations 
also involve the electric field which is determined by Eq. (11), and W° which is obtained 
from E using Eq. (9.b). For W, Eq. (12) would imply W oc S(k°), meaning that W is 
time independent. But Eq. (12) is only an approximation valid when k° 3> g 2 T. Thus 
from Eq. (12) one can only conclude that W evolves more slowly than . We will 
now see that the characteristic frequency of W is of order g 2 T. 

The appearance of the frequency scale g 2 T, which so far has not been discussed in 
the literature, is immediately obvious if we choose the z-axis in Eq. (9) parallel to k 
and expand W(v) in spherical harmonics Y im (v). Then the factor v • k in (9.c) equals 
k° — |k|cos6* v and does not mix W\ m with different m. Furthermore, the collision 
term is diagonal (cf. Eq. (5)), so that the W\ m with \m\ > 2 are completely decoupled 
from the gauge fields. The dynamics of these modes is thus governed by Eq. (9.c) 
but without the electric field. They perform oscillations with frequencies determined 
by v ■ k ~ g 2 T, which are driven by the noise term £, and which are damped by the 
collision term at a rate of order C ~ g 2 T. 

To include the \m\ < 1 modes of W and the gauge fields in this picture it is convenient 
to solve Eq. (9.c) for W, 



W(v) = I G k (v, v') \V ■ E + £(v')l . (14) 
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Here Gkiy, v') denotes the v-space inverse of the operator C — iv ■ k. It has cuts in the 
lower half of the complex /c°-plane which reflects the damping caused by the collision 
term. To understand the role of the electric field in Eq. (14) we insert this result into 



Since we consider ko ~ |k| ~ g 2 T, both terms on the lhs of Eq. (15) are of order g 2 T~E, 
while m^GkE on the rhs is of order TE. Consequently the lhs of Eq. (15) can be 
neglected completely when ko ~ g 2 T. This is a crucial point which will later allow us 
to drop the term DqE for all frequencies smaller than gT, even though both terms on 
the lhs of Eq. (15) are of the same order of magnitude when k ~ g 2 T. It implies that 
the two terms on the rhs of Eq. (15) must cancel. The electric field is thus entirely 
determined by the noise £ in such a way that the current m^W vanishes at leading 
order in g. Gauss' law (9.b) , together with Eq. (14) gives W° ~ g 2 W. Therefore both 
W° and W are small compared to W when ko ~ g 2 T. 

Finally we briefly recall the case k <C g 2 T which has been studied in great detail 
[6], [21]-[23], [16]. Again, Gauss' law implies that one can neglect W° in Eq. (9.c). 
Furthermore, one can approximate v ■ k ~ — v • k. This corresponds to dropping D W 
in Eq. (l.b). Then W is not dynamical, but it is fixed by the gauge fields and the 
noise at the same instant of time. For ko <C g 2 T the "magnetic" term /c^k x k x E in 
Eq. (15), which so far has not played any role, becomes relevant. It is now much larger 
than the kinetic term /c E which can be neglected. Then Eq. (15) gives ko ~ g A T as 
the characteristic frequency of the magnetic sector. 

To summarize the discussion of the linearized equations of motion (9), we have found 
that the characteristic frequency of the electric field E and the 4-current is given 
by uj p \ ~ gT. The characteristic frequency of W, i.e., of the I > 2 components of W is 
k° ~ g 2 T. Finally, the characteristic frequency of the magnetic fields is of order g A T. 

Now consider the effect of interactions in Eq. (1). Since D ~ V it is clear that none 
of the above order of magnitude estimates is changed by replacing V — > D. The basic 
picture of plasmon oscillations and multipole oscillations is unaffected, except that they 
occur in a quasi-static gauge field background. 

The effect of replacing d by D is less obvious and it depends on which frequencies 
are involved. First consider the plasmon oscillations. One can estimate the size of 
A using VA° ~ E. Because of equipartition E(x) and B(x) are of the same order 
of magnitude. This gives the same estimate for A°(x) as for A(x), i.e., A°(x) ~ gT. 
Consequently one can neglect gA° in the covariant time derivative acting on fields with 
frequencies of order gT, and Eq. (13) and the corresponding result for E are also valid 
in the interacting theory. Since Gauss' law contains a spatial derivative, the result for 
W° will look different in the presence of interactions but again the order of magnitude 
of W° is unchanged. 

Now consider A; < g 2 T, ignoring for a moment the physics of ko ~ gT. To estimate 
A we need to know the size of E(fc). We have seen that for A; <g 2 T the electric field 
is of the same order of magnitude as £. The latter can be estimated from the Fourier 



Eq. (10), 




(15) 
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transform of Eq. (4), which gives £(k) ~ A; ~ 1/2 |k|- 3 / 2 . Then one can use VA°(x) ~ 
/ d 4 ke- ik - x E(k), which gives A°(x) ~ g 2 T, g 3 T for k° ~ g 2 T, g 4 T, respectively. Note 
that these estimates are smaller than the one obtained from the equipartition argument 
above. This is because the electric field at a given time is dominated by Fourier 
components with frequencies of order gT. Nevertheless, Aq is big enough to be able 
to change the estimate d ~ g 4 T to D ~ g 3 T when D acts on a field with frequency 
of order g 4 T. Still, this would be a factor g smaller than D ~ g 2 T, and neglecting D 
would still be justified. 

We are now in the position to tentatively write down an effective theory for the 
physics associated with fc < g 2 T which corresponds to integrating out E and W 1 as 
dynamical degrees of freedom. We have seen that for both k° ~ g 2 T and k° ~ g 4 T one 
can neglect the term DqE. The spatial components of Eq. (l.a) can thus be replaced 
by 

D x B = m 2 D W. (16.a) 

In Eq. (l.b) we were able to neglect W° both for k° ~ g 2 T and k° ~ g 4 T. Therefore 
one can drop W° altogether. We have also seen that for A; ~ g 2 T we could neglect 
W. For k° ~ g 4 T we were able to neglect all time derivatives on the lhs of (l.b), and 
in particular the term Z)oW. Thus for both k° ~ g 2 T and k° ~ g 4 T the term DqW 
can be neglected and we can replace Eq. (l.b) by 

3(ci + v • D)v • W + (C + v ■ D)W = v • E + f , (16.b) 

where C\ is the I = 1 eigenvalue of the collision operator C. Eq. (16. a) is no longer a 
dynamical equation. Instead, it fixes the 3-current in terms of the gauge fields at the 
same instant of time. 

The only remaining question is whether plasmon oscillations affect the low frequency 
(k° <g 2 T) dynamics through interactions. To address this issue we consider the gauge 
field polarization tensor U. ftv (k) in the theory (1) at one loop. We are interested in 
ko < g 2 T, and loop momenta with q° ~ gT. Inside the loop one can neglect the effects 
of C and £ since q° ^> g 2 T. Without these terms Eq. (1) has the same form as the 
non-abelian Vlasov equations [5] which describe the hard thermal loop effective theory. 
Therefore the calculation of n Miy (/c) is precisely the same as the one in Ref. [11] where it 
was found that the leading order in g contribution is due to space-like loop momenta. 
This shows that the time-like loop momenta q° ~ gT, |q| ~ g 2 T considered here do 
not contribute at leading order in g. 

Written in A = gauge 4 , Eq. (16) is a Langevin equation which is purely dissipative, 
i.e., it contains only first order time derivatives of the dynamical degrees of freedom A 
and W, 

^ (v • A + W) = -(C + v • D) (3v • W + W) + f . (17) 

(J L 

4 Note that we have not assumed this gauge when deriving Eq. (16). 
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For the subsequent discussion I remind the reader of some features of a Langevin 
equation for some set of degrees of freedom tp a [25], 

J^a(*) = "/«[¥>(*)] +£«(*), (18) 

where / is a functional of tp at time t. In the present context the index a would 
represent the spatial coordinates, v and vector as well as color indices. Both Eq. (6) 
in A = gauge, and Eq. (17) are of this form, as well as the Langevin equation of 
Ref. [16]. In our case / is a local functional of tp, that is, it contains only the tp a and a 
finite number of spatial derivatives. In contrast, in the Langevin equation of Ref. [16] 
/ is a spatially non-local functional of the gauge fields 5 . Solving Eq. (18) for any 
given realization of the noise £ one generates an ensemble of field configurations for 
any time t. If the noise is Gaussian and white (that is, it has a frequency independent 
spectrum), 

(UtW))= 2 ™apS(t-t'), (19) 

the probability distribution for field configurations P{tp,t) satisfies the Fokker-Planck 
equation 



d 

dtp a 



d 

fa + TVt a p 



dp? 



P = 0. (20) 



An important case is that / is related to the derivative of a Hamiltonian H through 

riff 

fa = n a p — . (21) 

d<Pp 

Then the thermal equilibrium distribution P cq = exp[— H/T] is a stationary solution 
of the Fokker-Planck equation. Furthermore, for any initial configuration P(tp,t) ap- 
proaches P cq for large times. 

For Eq. (6) we have VL a p = 7~ 1 <5 Q , / 3 and f a = Vt^^dH^/dtpp) with the Hamiltonian 



= — 

m 2 



J d 3 xB 2 . (22) 



Thus for Eq. (6) the probability distribution approaches exp[— H 3d /T] for large times. 
Note that H 3d /T is the action of magnetostatic Yang-Mills theory [26] which is obtained 
for equilibrium quantities by dimensional reduction and by integrating out the Aq field. 

Beyond the leading and next-to- leading log approximation, for which Eq. (1) reduces 
to Eq. (6), it has so far not been understood whether Eq. (1) reproduces the correct 
thermodynamics of long distance Yang-Mills fields. The noise correlator of Ref. [16] 
depends on the gauge fields, and is therefore not of the form (19). With a field depen- 
dent noise correlator the time discretisation of the Langevin equation is ambiguous. 

5 In order to avoid confusion I would like to stress the following. By reintroducing the VF-field the 
Langevin equation of Ref. [16] can be written in a local form. Then, however, the equation for W 
would be a constraint and not an equation of motion. 
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In Ref. [16] this ambiguity was fixed by postulating that the Langevin equation yields 
the equilibrium distribution exp[— H 3A /T]. The considerations of Ref. [27] did not take 
into account the non-linear character of Eq. (l.b). 

We will now determine the equilibrium distribution which is generated by Eq. (17). 
Obviously the rhs of Eq. (17) can not be written like in Eq. (18) with an / of the form 
(21). It is well known, however, that the equilibrium distribution does not specify the 
function f a uniquely [25]. One may add an extra term F a , 

OH 

f a = VL a p- h F a , (23) 

without changing the equilibrium distribution provided that F a satisfies 

d(f a dif a 

We will see that our Langevin equation (17) is indeed of the form (18), (23) with an 
F satisfying the condition (24). From Eq. (4) we see that we have to identify 

Q Q/3 <-> \c{\, v')5(x - *!)5 ab . (25) 

The next question concerns the relevant Hamiltonian. Eq. (1) was obtained from the 
Hard Thermal Loop effective theory for which the Hamiltonian reads [28] 

#htl = \ J d'x {E 2 + B 2 + ml jf H/ 2 } . (26) 

The effective theory (16) is obtained from Eq. (1) by integrating out E and W^. Thus 
one can expect that the relevant Hamiltonian is obtained from (26) by dropping these 
fields, i.e., 

H = l -Jo?x^B 2 + mlJ^W 2 Y (27) 

Using 

- = DxB, W = <W, (28) 

together with Eq. (25) we indeed obtain the terms on the rhs of Eq. (17) which contain 
the collision term. The remaining terms on the rhs of Eq. (17) have to be identified 
with F a , 

F a —v ■ D (3v • W + W^j . (29) 

We will now see that for this F Q both the lhs and the rhs of Eq. (24) are zero, so that 
Eq. (24) is indeed satisfied. First consider the lhs of Eq. (24). When the (^-derivative 
acts on the gauge field contained in v-D one obtains zero due to the contraction of color 
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indices and the antisymmetry of the structure constants. When it acts on 3v ■ W + W 
the result is an integral over a total spatial derivative which again vanishes. The rhs 
of Eq. (24) can be written as 



with ip = v • A + W. This again is an integral of a total derivative and vanishes. 

We conclude that the probability distribution generated by the Langevin equation 
(17), with W given by Eq. (16. a), approaches for large times the Boltzmann distribution 
P eq = exp(— if/T) with the Hamiltonian (27). This is not the equilibrium distribution 
corresponding to dimensional reduction because it contains the additional field W. In 
H, however, W is not coupled to the gauge fields. Equal time correlation functions 
of A computed with the help of H are thus the same as in 3-dimensional Yang-Mills 
theory. 

Another interesting observation is the following. Without a kinetic term DqE there 
are no propagating gauge field waves in the theory (16). These would cause the same 
non-local UV divergences as in classical thermal Yang-Mills theory [20] , which are still 
present in Eq. (1) since the effect of the W fields can be neglected in the ultraviolet. 
Therefore the UV divergences of Eq. (16) can be expected to be local. 

To summarize, we have found that Eq. (1) describes plasmon oscillations, multipole 
fluctuations of color charge, and the non-perturbative gauge field dynamics. These 
processes are associated with characteristic frequencies gT, g 2 T, and g 4 T, respectively. 
An effective theory (16) was constructed which reproduces the slow (ko ~^g 2 T) dynam- 
ics at leading order in g. Previously [6, 22, 16, 12] it was implicitly assumed that all 
modes with k ^> g A T decouple from the non-perturbative gauge field dynamics. Here 
it was shown that this is indeed the case for the plasmon oscillation. To see whether 
the same is true for the multipole fluctuations requires a more detailed analysis. 

Acknowledgements I would like to thank P. H. Damgaard, M. Laine, G.D. Moore 
and K. Rummukainen for useful discussions and suggestions. This work was supported 
by the TMR network "Finite temperature phase transitions in particle physics", EU 



contract ERBFMRXCT97-0122. 

References 

[1] A. D. Linde, Phys. Lett. B96 (1980) 289; D. J. Gross, R. D. Pisarski, L. G. Yaffe, 
Rev. Mod. Phys. 53 (1981) 43. 

[2] For reviews of high temperature dimensional reduction see, e.g., M. E. Shaposh- 
nikov, hep-ph/96 10247; A. Nieto, Int. J. Mod. Phys. A12, 1431 (1997) [hep- 
ph/9612291]; M. Laine, hep-ph/9707415. 

[3] For a review of electroweak baryon number violation and baryogenesis see, 
e.g., V.A. Rubakov, M.E. Shaposhnikov, Usp. Fis. Nauk 166 (1996) 493 [hep- 
ph/9603208]. 




9 



[4] E. Braaten and R. Pisarski, Nucl. Phys. B337 (1990) 569; Phys. Rev. D45 (1992) 
1827; J. Frenkel and J.C. Taylor, Nucl. Phys. B 334 (1990) 199; J.C. Taylor and 
S.M.H. Wong, Nucl. Phys. B346 (1990) 115. 

[5] J. P. Blaizot, E. Iancu, Nucl. Phys. B417, 608 (1994) [hep-ph/9306294]; V.P Nair, 
Phys. Rev. D48, 3432 (1993) [hep-ph/9307326]. 

[6] D. Bddeker, Phys. Lett. B426, 351 (1998) [hep-ph/9801430]. 

[7] P. Arnold, D. T. Son, L. G. Yaffe, Phys. Rev. D59, 105020 (1999) [hep- 
ph/9810216]. 

[8] M. A. Valle Basagoiti, hep-ph/9903462; 

[9] D. F. Litim and C. Manuel, Phys. Rev. Lett. 82, 4981 (1999) [hep-ph/9902430]; 

[10] J. Blaizot and E. Iancu, Nucl. Phys. B557, 183 (1999) [hep-ph/9903389]. 

[11] D. Bodeker, Nucl. Phys. B566, 402 (2000) [hep-ph/9903478]. 

[12] P. Arnold, L. G. Yaffe, Phys. Rev. D62, 125014 (2000) [hep-ph/99 12306]; Phys. 
Rev. D62, 125013 (2000) [hep-ph/9912305]. 

[13] The lhs of Eq. (6) can be estimated as (g 2 T) 2 A, while on the rhs 7E ~ 
T/ \og(l/g)k A. This gives k ~ log(l/^ 4 T. 

[14] G. D. Moore, Nucl. Phys. B568, 367 (2000) [hep-ph/9810313]. 

[15] G.D. Moore, Phys. Rev. D62, 085011 (2000) [hep-ph/0001216]. 

[16] P. Arnold, Phys. Rev. D62, 036003 (2000) [hep-ph/9912307]. 

[17] G.D. Moore, C. Hu, B. Muller, Phys. Rev. D58, 045001 (1998) [hep-ph/9710436]. 

[18] G. D. Moore, K. Rummukainen, Phys. Rev. D61, 105008 (2000) [hep-ph/9906259]. 

[19] D. Bodeker, G.D. Moore, K. Rummukainen, Phys. Rev. D61, 056003 (2000) [hep- 
ph/9907545]. 

[20] D. Bodeker, L. McLerran, A. Smilga, Phys. Rev. D52, 4675 (1995) [hep- 
th/9504123]. 

[21] P. Arnold, D. Son, L. G. Yaffe, Phys. Rev. D55, 6264 (1997) [hep-ph/9609481]; 
P. Huet, D. T. Son, Phys. Lett. B393, 94 (1997) [hep-ph/9610259]. 

[22] D. BSdeker, Nucl. Phys. B559, 502 (1999) [hep-ph/9905239]. 

[23] P. Arnold, D.T. Son, L.G. Yaffe, Phys. Rev. D60, 025007 (1999) [hep- 
ph/9901304]. 

[24] E.M. Lifshitz, L.P. Pitaevskii, Physical Kinetics (Pergamon Press, Oxford 1981). 



10 



[25] For an introduction and review, see J. Zinn- Justin, Quantum Field Theory and 
Critical Phenomena, Clarendon Press, Oxford, 1993. 

[26] E. Braaten and A. Nieto, Phys. Rev. D 53 (1996) 3421 [hep-ph/95 10408]. 

[27] D. F. Litim and C. Manuel, Phys. Rev. D 61, 125004 (2000) [hep-ph/99 10348]. 

[28] J. Blaizot and E. Iancu, Nucl. Phys. B421, 565 (1994) [hep-ph/9401211]; 
V.P Nair, Phys. Rev. D50, 4201 (1994) [hep-th/9403146]. 



11 



